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The study of core noise from turbofan engines has become more important as noise from other sources such as the 
fan and jet were reduced. A multiple-microphone and acoustic-source modeling method to separate correlated and 
uncorrelated sources is discussed. The auto- and cross spectra in the frequency range below 1000 Hz are fitted with a 
noise propagation model based on a source couplet consisting of a single incoherent monopole source with a single 
coherent monopole source or a source triplet consisting of a single incoherent monopole source with two coherent 
monopole point sources. Examples are presented using data from a Pratt & Whitney PW4098 turbofan engine. The 
method separates the low-frequency jet noise from the core noise at the nozzle exit. It is shown that at low power 
settings, the core noise is a major contributor to the noise. Even at higher power settings, it can be more important 
than jet noise. However, at low frequencies, uncorrelated broadband noise and jet noise become the important 
factors as the engine power setting is increased. 
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Nomenclature 

amplitude of a coherent monopole point 
source A 

amplitude of a coherent monopole point 
source B 

resolution bandwidth, B e = l/T d = 
r s /NP = 11.71875 Hz 
speed of sound, m/s 

function providing solution space boundary 
on SPLjy 
expected value 

two-sided Fourier transform operator 
frequency 

upper frequency limit, f c = 1/ (2A t) = rJ2, 
24,000 Hz 

cross power spectral density function 
defined for nonnegative frequencies only 
(one-sided) 

auto power spectral density function 

defined for nonnegative frequencies only 

(one-sided) 

imaginary unit, </--l 

wave number, 2nf /c D , m 1 

number of frequencies, f c / A / = N/2 (2048) 

molecular weight of air, 28.97, kg/k mol 

amplitude of uncorrelated monopole point 

noise source, N 

segment length, number of data points per 
segment (4096) 

number of disjoint (independent) segments, 
lid Bet{ otal 234 

number of overlapped segments/blocks 
probability confidence interval, % 
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pressure wave function 
20 juPa 

distance from the point source, m 
gas constant, 83 1 4 J/k mol K 
distance from the point source a to 
microphones i, m 
distance from the point source a to 
microphones j, m 
autocorrelation function 
cross-correlation function 
radial distance, m 

sample rate, 48,000 samples per second 
0.3045 m/ft 

two-sided autospectral density function 
two-sided cross-spectral density function 
source strength of sources A, B, and N 

autospectra error cost function 
cross-spectra magnitude error cost 
function 

cross-spectra phase-angle-error cost 
function 

total cost function 
record length of the segment 
total record length, »20 s 
air temperature, 80° F 
temperature, K 
time, s 

functions providing solution space 
boundaries on x A and x B 
x position of sources A and B 
spatial coordinates of microphones and 
sources 

high and low boundaries on x A 
high and low boundaries on x B 
frequency resolution, 1 /T d , 11.718 Hz 
sampling interval, 1 / r s , 1/48,000 s 
Dirac delta function 
ratio of specific heats, 1.4 
magnitude-squared coherence function of 
noise 

magnitude-squared coherence 
function 

time displacement, s 
phase angle, rad 
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Subscripts 

o 

= some arbitrary specific value 

X 

= signal x 

y 

= signal y 


I. Introduction 


U NDERSTANDING turbofan engine noise is a key priority of 
NASA. Consequently, diagnostic procedures to identify 
dominant sources and changes in dominant sources are being 
developed. This paper discusses a diagnostic procedure that 
separates correlated far-held turbofan noise from uncorrelated far- 
held turbofan noise using a multiple-microphone method and point- 
source acoustic models. The auto- and cross spectra in the frequency 
range below 1 000 Hz are htted with a noise propagation model based 
on a source couplet consisting of a single incoherent source with a 
single coherent source or a source triplet consisting of a single 
incoherent source with two coherent point sources. Two noise 
signals are coherent if they can be aligned with one another so that the 
coherence calculated using a periodogram averaging method is 
greater than the coherence of two random signals using the same 
periodogram averaging method [JJ. As fan noise and jet noise from 
turbofan engines are reduced, the issue of core-noise reduction 
becomes more important. Although core noise may be reduced by 
acoustic liners, the issue of measuring the amount of reduction 
becomes increasingly important. The proposed scheme separates 
coherent noise from random jet noise and should enable the 
measurement of the effectiveness of core-noise reduction liners. 
The scheme was developed as part of a research program to study 
core noise from a Pratt & Whitney PW4098 turbofan engine 
(Miles [1-3].) 

Procedures using coherence-based techniques have been 
developed for extracting acoustic signals buried in noise. The 
coherent output power method for noise-source identification is 
discussed by Bendat and Piersol [4—6] . The application of this 
technique that is of interest is the use of coherent output power 
spectra to separate and identify correlated combustion noise in far- 
field measurements of turbofan engine noise. Karchmer [7] and 
Karchmer et al. [8], use the coherence function calculated from 
internal microphone measurements of fluctuating pressures in the 
combustor and far-field acoustic pressures to determine the 
correlated combustion noise of a YF102 turbofan engine at far-field 
locations by calculating the coherent output power spectrum. 

The three-signal coherence technique was developed by Chung 
[9, HI] for flow noise rejection. A similar technique was developed 
and used by Krejsa [11], The three-signal coherence technique was 
used by Shivashankara [17] to study core noise in a Pratt & Whitney 
JT9D. It was used by Hsu and Ahuja [ 13 ] to separate ejector internal 
mixing noise from far-field measurements and by Stoker et al. [_14] to 
separate wind-tunnel background noise and wind noise from 
automobile interior noise measurements. It was used by Michalke 
et al. [ 15 ] to study sound in a circular duct with mean flow. The 
method is also discussed by Minami and Ahuja [16]. 

The inherent coherent properties of jet and core noise have been 
used in several source-separation procedures. The radiated field of jet 
noise has a low coherence when measured at two widely separated 
points. Parthasarathy et al. [17] attributed this to the fact that the jet- 
noise sources are in motion so that the source frequencies undergo 
large Doppler shifts as the noise is radiated to the far-field. Core noise 
has a high coherence when measured at two widely separated points. 
This is attributed to the fact that the frequencies of the radiated core 
noise are preserved unchanged in the far field (Parthasarathy et al. 
[17]). Consequently, three-signal coherence technique is especially 
useful to study turbofan core noise using widely spaced 
microphones. 

A method of separating jet noise and core noise using 
autocorrelations and cross correlations was developed by 
Parthasarathy et al. [17]. The model presented assumes that the 


microphones are located on the arc of a circle about a single source 
having a known position. This is a couplet source, which produces a 
core-noise sound radiation field that is correlated over wide 
microphone spacing and a jet-noise sound radiation field that is not 
correlated over wide microphone spacing. Because the microphones 
are on the same arc, the model cross correlations calculated have zero 
propagation time delay. In addition, a model for the jet-noise 
directivity is used. The resulting system of three unknowns and three 
equations can be solved exactly. In applying this method to data from 
a small jet, the source location is known and a correction to the 
amplitude correlation measurements is made so that the microphones 
are located on the arc of a circle about a source having a known 
position. In applying this method to measured data, I assume the 
cross correlations used by Parthasarathy et al. were also time-shifted 
to remove the propagation time delay. An extension of the method to 
obtain spectral information is also discussed using nonzero values of 
propagation time delay. The phase-angle shift due to propagation 
time delay is not discussed. Again, this seems to indicate that the 
propagation time delay is being set to zero. The method calculates the 
jet-noise spectrum, the core-noise spectrum, and the directivity ratio 
of the core noise at a particular angle as a function of frequency from 
measured auto- and cross spectra. The method was applied to 
experimental data obtained in a small-scale facility. The spectral 
formulation of this method was used by Tesson [J_8, HI] to study jet 
noise and core noise from the static test of a small gas turbine engine 
in an anechoic chamber. 

A method to achieve separation of two different correlated noise 
sources from far-field measurements contaminated by extraneous jet 
noise using multiple microphones was developed by Minami and 
Ahuja [16], The equations discussed use measured autospectra and 
cross spectra. The model assumes that the source noise can be 
represented by a triplet source consisting of a correlated noise source 
A, correlated noise source B, and an uncorrelated noise source. In this 
analysis by Minami and Ahuja, it is assumed that all the sources are at 
a single triplet location; however, in the formulation presented, no 
source location or microphone location information is used. The five- 
microphone method described involves solving a set of 55 equations 
for 55 unknowns at each frequency of interest. At each microphone 
location, the spectrum of correlated noise source A, correlated noise 
source B, and an uncorrelated noise source is obtained as a function 
of frequency. The method was validated using analytical 
simulations. 

Beamforming and phased-array methods used in acoustic testing 
require similar measurements and use a similar source description 
[20-26]. However, they use more than four microphones and are 
based on a fixed geometry, whereas the procedure described herein 
uses a floating geometry. A phased-array procedure using a linear 
array of 1 4 microphones was used to study low-frequency noise from 
a Rolls-Royce/SNECMA Olympus engine by Flynn and Kinns [ 27 ] 
and Billingsley and Kinns [281. Concurrently, the polar correlation 
technique was developed by Fisher et al. [29] and used to study 
Rolls-Royce Viper 60 1 engine noise and Rolls-Royce RB2 1 1 engine 
noise. 

Source characterization by correlation techniques used by Blake 
and Lynch [ 30 ] to study trailing-edge noise and fan noise require 
similar measurements. A procedure for source characterization using 
a closely spaced pair of microphones in the far field was developed by 
Kinns [31], 

As pointed out by Williams [32] in a discussion of antisound, there 
is no unique specification for the origin of a given sound field. One 
may only think in terms of an equivalent source region that shares the 
measurable far-field properties of the real sources. As he states, 
“There is no difference between the sound of a naturally occurring 
phase-related multiple-source system and that created artificially 
with control (electronic and acoustic apparatus) but the latter can be 
changed and adapted at the ‘flick of a switch' as it were.” The source- 
separation method described herein uses this concept. The particular 
source models described were chosen for convenience and have a 
feasible analogy. The chosen mathematical models permit the 
creation of the correct directivity and phase relationships between the 
microphone signals as a function of frequency. 
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The method presented herein also models the auto- and cross- 
spectrum measurements made with multiple microphones. Rather 
than use a single mathematical-equation-based model, the method 
uses a sequence of models. Each model consists of a collection of one 
or more point-source solutions to a three-dimensional acoustic wave 
equation in spherical coordinate form. Each element of the sequence 
of models is based on propagation from a collection of point-source 
sound radiators on the engine axis that can have arbitrary positions on 
the engine axis. Consequently, the microphones need not be on a 
circular arc at the same distance from the source. In addition, for each 
source, the directivity is the same in each direction for each source. 
The method tests the measured data at each frequency against various 
models. The first model in the sequence assumes an incoherent 
source of strength SPL N (f) at the location x N = 0.0 and y N = 0.0. 
The second model in the sequence assumes that the source is a 
doublet consisting of an incoherent source of strength SPL N (f) at the 
location x N = 0.0 and y N = 0.0 and an additional coherent source of 
strength SPL A (/) at location x A (f) with v A = 0.0. The third model in 
the sequence assumes that the data are the result of a source triplet 
consisting of an incoherent source of strength SPL N (f) at the 
location x N = 0.0 and y N = 0.0, a coherent source at location x A (f) 
with y A = 0.0 of strength SPL A (/), and an additional coherent 
source at location x B (f) with y B = 0.0 of strength SPL s (/). 

For each model in the sequence, the appropriate parameters were 
adjusted to give the best fit to the data at each frequency. The results 
presented here are those obtained by the method when one assumes 
that the engine noise source can be modeled with a source doublet or 
triplet. 

One part of a source doublet and triplet consists of an incoherent 
source of strength SPL W (/) at the location x N = 0.0 and y N = 0.0. 
The source doublet has an additional coherent source of strength 
SPL a (/) at location x A (f) with y A = 0.0. The source triplet has two 
additional coherent sources at location x A (f) with y A = 0.0 and 
x B (f) with y B = 0.0. As part of the procedure, a comparison of 
results obtained with the use of a source doublet with the results 
obtained using a source triplet is made. A source position is assumed 
to have a y coordinate of zero and the x coordinate along the engine 
axis is found as a function of frequency. This is done because low- 
frequency noise can be radiated to the far field from the inlet, the 
engine case, the exhaust nozzle, and the jet. To solve the resulting 
acoustic equations for a small set of parameters, a solution method 
that provides an optimum solution in a least-squares sense without 
derivatives was used. The method is set up to do a least-squares curve 
fit to match the autospectra and the cross-spectra magnitude and 
phase measurements between all microphones. The method assumes 
the existence of one or more point sources that can produce the same 
measurements. The problem then becomes one of finding the 
appropriate point sources. The current strategy is to have the curve fit 
adjust the amplitude and x location of the sources that determine the 
cross spectrum at each microphone and the amplitude of a random- 
noise source. For the single coherent source case, three parameters 
are found at each frequency. When a model with two coherent 
sources is used, one has five parameters to determine at each 
frequency. Examples are presented calculated with a four- 
microphone array using data from a Pratt & Whitney PW4098 
turbofan engine. 

II. Analysis Method 

The core noise is assumed to be propagating in acoustic modes in 
the turbofan engine. The coherent acoustic energy leaves the nozzle 
and travels through the turbofan engine shear layer to a ground 
microphone. Additional acoustic energy from the jet and from 
random scattering reaches the same microphone. In this test program, 
four microphones at 150 ft and angular position of 100, 110, 120, and 
130 deg measured from the inlet were used. Consequently, the 
measurements available at each test condition are four sound- 
pressure spectrums and six sound-pressure cross spectrums, each 
consisting of a magnitude and a phase angle. These 1 6 measurements 
are available as a function of frequency. Four real variable and six 
complex variable acoustic equations relate the measurements and the 


model parameters. The angular separation of at least 10 deg means 
that the jet noise measured at any microphone can be assumed to be 
incoherent above some frequency with the jet noise at any another 
microphone. 

The basic procedure is to assume that some combination of 
coherent sources and an incoherent source will produce the measured 
auto- and cross spectra. Coherent sources are assumed to be on the 
turbofan engine centerline at some x(f ) position and with v = 0. The 
incoherent source is assumed to be at x = 0 and y = 0. Using the 
results in Appendix A, in which the one-sided spectrum functions for 
acoustic signals from a point source are derived, the acoustic model 
equations shown in Appendices B, C, D, and E are defined. The 
following models were considered: 

1) The two-parameter model has a coherent source of magnitude A 
at x = x a and y = 0 (Appendix B). 

2) The three-parameter model has a coherent source of magnitude 
A at x = x a and y = 0 (Appendix C) and an incoherent source of 
magnitude TV at x = 0 and y = 0. 

3 ) The four-parameter model has a coherent source of magnitude A 
at x = x a and y = 0 and a coherent source of magnitude B at x = x b 
and y = 0 (Appendix D ). 

4) The five-parameter model has a coherent source of magnitude A 
at x = x a and y = 0 and a coherent source of magnitude B at x = x b 
and y = 0 with an incoherent source of magnitude TV at x = 0 and 
y = 0 (Appendix E). 

Consequently, the number of equations greatly exceeds the 
number of model parameters used in the study discussed herein. Note 
that the radiation pattern assumed for these models is one that 
radiates equally in all directions. These are monopole sources. 

The measured autospectra and cross spectra required as input to 
solve the acoustic model equations is experimentally determined and 
subject to nominal experimental error and statistical uncertainties. 
Also, the measured autospectrum includes random noise in addition 
to coherent signals from propagating waves. Consequently, a 
solution method that provides an optimum solution in a least-squares 
sense without derivatives was used. Algorithms for minimization of a 
cost function without derivatives are discussed by Brent [33], The 
search technique used in this study is described by Powell [34], and a 
Fortran computer code for this algorithm is given by Shapiro and 
Goldstein [ 35 ] and Kuester and Mize [36]. The code used was a 
modified version of the one in Shapiro and Goldstein [35], which was 
updated to be in a Fortran-90 style. 

A typical cost function is discussed in Appendix F. The cost 
function serves three purposes. In addition to providing a function 
defining the solution, the cost function is also used to break the 
symmetry for the case in which the solution has two coherent 
sources. A plot of the locations x A and x B can be confusing for this 
case because the coherent monopole sources are symmetric. 
Consequently, at a particular frequency, interchanging positions and 
magnitudes does not change the answer. However, the plots of the 
locations x A and x B can look strange. To explicitly break the 
symmetry, a penalty function is used to force the condition x B < x A . 

The third function of the cost function is to confine the solution to a 
reasonable domain. After looking at the results of many 
unconstrained searches, two search regions were selected based on 
frequency. For /<164Hz, we use ll<x A <100 and 
—25 < x B < 100. However, for / > 164 Hz, we use 1 1 < x A < 25 
and —25 < x B < 11 (where the locations are in units of feet). The 
resulting penalty functions prevent convergence to a solution in 
which x A = 3000. 

III. Experimental Setup 

To demonstrate the usefulness of the procedure for separating 
correlated and uncorrelated noise sources, measurements made in the 
far field of a Pratt & Whitney PW4098 turbofan engine will be used. 
The measurements were made in a study of aircraft turbofan engine 
core noise conducted as part of the NASA Engine Validation of 
Noise Reduction Concepts (EVNRC) program. 

The spectral estimate parameters are shown in Table 1 . The signal 
processing algorithms used were written in Fortran. They are based 


64 


MILES 


Table 1 Spectral estimate parameters 


Parameter 

Value 

Segment length (data points per segment), NP 

4096 

Sample rate r, samples per second 

48,000 

Segment length, T d = NP/ r, s 

0.08533 

Sampling interval, At = \/r , s 

2.0833 x nr 5 

Frequency step, A / = \/T d , Hz 

11.718 

Upper frequency limit, f c = 1 /2A t = r/2, Hz 

24,000 

Number of frequencies, Ly = f c / A / = N p/2 

2048 

Time delay, r = 6323/48, 000 s 

0.1317 

Number of independent samples 

234 

Overlap 

0.50 

Sample length, s 

20 


on subprogram modules developed by Steams and David [37]. In the 
calculations, the segments were overlapped by 50%. Figure J_ shows 
the angular placement of the far-field microphones on a 150-ft radius. 
This analysis uses the microphones at 100, 110, 120, and 130 deg 
measured from the inlet. The coordinate system used herein has the x 
axis along the engine centerline, increasing to the right. The y axis is 
perpendicular and increases toward the top of the plot. The center 
(x = 0 and y = 0) is at the engine center. 

In calculating the cross spectra between the microphones, no time 
delay is used. Consequently, the cross spectra between the 
microphones have a phase shift related to the time a wave front from a 
particular source reaches each microphone. 


independent segments/blocks n d , we take n d = n a , where n a is the 
number of overlapped segments/blocks. The value of n a used herein 
is n a = 460. Consequently, fromEq. (3), we have y% n = 0.0065. The 
coherence threshold y„ n is discussed by Carter [38,39], Halliday et al. 
[40], and Brillinger [41_] . The coherence threshold y^ n has a value that 
is greater than 95% of the values of the coherence of two independent 
time series calculated using n d disjoint data segments/blocks. The 
coherence threshold used herein is calculated using the number of 
overlapped segments, n a . As part of the set of figures for each case 
used to show the results, the coherence and the threshold coherence 
yj[ n will also be shown. The data used herein had coherence values 
above the threshold coherence y^ n . 


B. Cross-Spectra Phase-Angle Sampling Errors 

In Bendat and Piersol [6] and Piersol [42], the random error in the 
phase estimates due to statistical sampling is given in terms of the 
standard deviation of the estimated phase angle 6 n by 


0012 (/)] 


sin 1 


([i-Awn 

| \YnW2K j 


(4) 


where a[9 n (f)\ is measured in radians, and as used herein, n a is 
selected to be the number of overlapped segments or blocks used in 
the spectral calculations. For the special case in which the term in 
braces is small, Eq. (4) becomes 


IV. Results 

A. Cross-Spectra Validity 

The signals from each pair of microphones used to calculate the 
cross spectra only produce valid measurements if the coherence is 
greater than some threshold coherence. The coherence function is 
given by 


rieM) 


|Gfl,9 2 (/)| 2 


(1) 


The measured coherence calculated using segments overlapped by 
50% is given by 


Yle,<f) = 


iE£t^o-./)^o-./)r 
E"=i l^,(i,/)| 2 E"=r \Xe 2 (i,f)\ 2 


( 2 ) 


In Miles [j_], comparisons were made of a coherence threshold based 
on aligned and unaligned coherence and one based on analytical 
coherence threshold using computer simulation. Results were 
obtained from computer simulation that show good agreement with 
the theoretical estimate of the analytical coherence threshold: 

Ynn = 1 — (1 — (3 ) 


where we use herein P = 0.95 and instead of the number of 



Fig. 1 Acoustic arena and microphone array at Pratt & Whitney test 
stand Cll, West Palm Beach, Florida, for EVNRC phase-2 tests. 


cr[e i2 (f)\ 


[ i - k ? 2 (/)] 1/2 

\Yu\V2th, 


(5) 


where for the unknown coherence y 2 2 (/), the estimated coherence 
Yg^(f) from Eq. (2) is used. A plot of the standard deviation of the 
phase angle in degrees versus coherence is shown in Fig. 2. 

When the coherence is greater than 0.15, Fig. 3 shows that the 
standard deviation should be less than 5 deg. Straight-line fits were 
made to the phase-vs-frequency data. The standard deviation of the 
error was greater than one would calculate from the coherence (i.e., 
around 10 deg when one calculated 5 deg). Consequently, the phase- 
angle measurements might be showing propagation effects due to 
wind shear or temperature gradients or a change in source location 
with frequency. In addition, the error-vs-frequency plots seemed to 
be in phase for each pair. As a consequence, instead of a using a fixed 
source location in the model, the y coordinate of a source is fixed at 
zero and the x coordinate is free to vary. 



Coherence 

Fig. 2 Standard deviation of phase angle of G 34 based on and n a . 
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Frequency, Hz 


Frequency, Hz 


Frequency, Hz 


d) e) f) 

Fig. 3 Three-parameter models at 150 ft with lVl COIT = 582 rpm : a) autospectrum 120-deg microphone, b ) cross-spectrum magnitude between 120- and 
130-deg microphones, c) cross-spectrum phase angle between 120- and 130-deg microphones, d) coherence between 120- and 130-deg microphones, 

e) SPL A and SPI, v parameters, and f) x A location parameter. 


C. Point-Source Models 

The method discussed was tried out over a range of operating 
conditions. For each case, all four models were used. Four typical 
cases presented herein are shown in Figs. 3-6 for lVl colr values of 
582, 1622, 1900, and 2743 rpm. Only one acoustic model result is 
shown for each case. The selected case uses the fewest number of 
parameters to achieve a good result. For each case, all four 
autospectra and six cross spectra are calculated. Flowever, only the 
following items are shown herein: 


1 ) The measured and calculated autospectrum is shown at 120 deg. 

2) The measured and calculated cross-spectrum magnitude is 
shown between the 120- and 130-deg microphones. 

3) The measured and calculated cross-spectrum phase angle is 
shown between the 120- and 130-deg microphones. 

4) The measured coherence and analytic threshold coherence is 
shown between the 120- and 130-deg microphones. 

5) The sound-pressure amplitude curve-fit parameters are shown. 

6) The correlated source locations are shown. 




Fig. 4 Five-parameter models at 150 ft with V 1 ,„ rr = 1622 rpm : a) autospectrum 120-deg microphone, b ) cross-spectrum magnitude between 120- and 
130-deg microphones, c) cross-spectrum phase angle between 120- and 130-deg microphones, d) coherence between 120- and 130-deg microphones, 
e) SPL A and SPI, V parameters, and f) x A location parameter. 
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Frequency, Hz Frequency, Hz 



Frequency, Hz 


d) e) f) 

Fig. 5 Five-parameter models at 150 ft with iVl corr = 1900 rpm : a) autospeetrum 120-deg microphone, b) cross-spectrum magnitude between 120- and 
130-deg microphones, c) cross-spectrum phase angle between 120- and 130-deg microphones, d) coherence between 120- and 130-deg microphones, 

e) SPL A and SPL W parameters, and f) x A location parameter. 




Fig. 6 Three-parameter models at 150 ft with /V 1 ,.„ rr = 2743 rpm: a) autospectrum 120-deg microphone, b) cross-spectrum magnitude between 120- 
and 130-deg microphones, c) cross-spectrum phase angle between 120- and 130-deg microphones, d) coherence between 120- and 130-deg microphones, 
e) SPL a and SPI. V parameters, and f) x A location parameter. 


In the modeling scheme used herein, for each case, the jet noise 
was initially assumed to be uncorrelated between the microphones at 
all frequencies. The results show that this assumption may be true at 
frequencies greater than 164 Hz. Furthermore, it is assumed the basic 
radiation pattern of the correlated noise as specified by the cross- 
spectrum phase angle could be represented by a single coherent point 
source at some x location. In all cases, a better result is obtained by 
using an acoustic model that assumes that an uncorrelated noise 
source is present and letting the computer solution procedure 
determine if its value is significant. In some cases, a second coherent 
point source gives a slightly better result. Only the coherent sources 


create the cross spectra. The autospectra are due to the coherent and 
incoherent sources. 

The general formulation is given in Appendix A. The general 
models are shown in Appendices B. C, D, and E 

1. N l cor r = 582 rpm 

At the engine-idle N l corr = 582-rpm operating condition, the core 
noise is dominated by a low-frequency combustion instability tone at 
327 Hz, as shown in Fig. 3a, in which the 120-deg autospectrum is 
shown. Results are shown for the three-parameter doublet- source 
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model ( A , N, and x A ). The noise radiation pattern can be represented 
by a single point source, as shown in Figs. 3b and 3c, in which the 
cross spectrum between the 120- and 130-deg microphones are 
shown. This cross spectrum is shown because the coherence 
calculated using the 120- and 130-deg microphones is above 0.04, as 
shown in Fig. 3d, and the phase-angle error is less than 9 deg (using 
Fig. 2). The results shown are for the three-parameter couplet-source 
model (A, N, and x A ) described in Appendix C. As shown in Fig. 3a, 
at a few low-frequency points, the computer solution included an 
uncorrelated noise source. The coherence shown in Fig. 3d is above 
the coherence threshold y% n = 0.0065, shown as a dashed line in 
Fig. 3d. However, the coherence shows a lot of variation above 
400 Hz, which signifies a greater statistical deviation in the phase 
angle and more uncertainty in the location solutions. The values of 
SPL A and SPL W are shown in Fig. 3e. At most points, the 
uncorrelated noise level was too low to plot, as shown in Fig. 3e. As 
shown in Fig. 3f, to match the measured phase angle, the computer 
tends to find solutions in which x A = 10 ft at low frequencies and 
20 < x A < 40 at higher frequencies. 

2. N\ co „ = 1622 rpm 

At the N\ con = 1622-rpm operating condition, results are shown 
for the five-parameter triplet-source model (A, B, N, x A , and x B ) 
described in Appendix E. As shown in Fig. 4a, this model does 
a good job of representing the 120-deg microphone measured 
autospectrum. The noise radiation pattern is well-represented by 
two point sources, as shown in Figs. 4b and 4c, in which the 
magnitude and phase of the cross spectrum between the 120- and 
130-deg microphones are shown. Below 400 Hz, this model does a 
fine job. The coherence between the 120- and 130-deg microphones 
shown in Fig. 4d is above the coherence threshold yTm = 0.0065. 
Below 400 Hz, the coherence is greater than 0.10. Consequently, the 
cross-spectrum phase angles look good below 400 Hz. 

The source solution SPL values are shown Fig. 4e, and the source 
solution location values are shown in Fig. 4f. The variation of the 
source-A location solution below 1 64 Hz corresponds to the change 
in thelocation of a low-frequency noise source from near the end of 
the potential core of the jet toward the nozzle as higher-frequency 
solutions are found. Consequently, we associate SPL A with the jet 
source in the frequency range below 164 Hz. 

The variation of the source amplitude of the solution for source B 
shown in Fig. 4e and the variation of the far-field SPL fi autospectra 
shown in Fig. 4a have the same frequency dependence that the 
combustion modes have, in that they seem to be peaking in the 0- to 
200-Hz band and the 200- to 400-Hz band. This characteristic 
behavior was reported in a study of PW4098 annular combustor 
modes by Miles [3]. This suggests that the correlated noise source B 
is associated with the combustor source in the 0- to 400-Hz frequency 
range. The source locations of the two point sources shown in Fig. 4f 
have an interesting symmetry that will be discussed in the next 
section. 

3. N l co „ = 1900 rpm 

Again, at the AH corr = 1900-rpm operating condition, results are 
shown for the five-parameter triplet-source model (A, B, N, x A , and 
x B ) described in Appendix E. As shown in Fig. 5a, this model does a 
good job of representing the 120-deg microphone measured 
autospectrum. 

The noise radiation pattern is well-represented by two point 
sources, as shown in Figs. 5b and 5c, in which the magnitude and 
phase of the cross spectrum between the 120- and 130-deg 
microphones are shown. The coherence shown in Fig. 5d is above the 
coherence threshold yl„ = 0.0065 and greater than 0.04 at most 
frequencies less than 500 Hz. Consequently, the cross-spectrum 
phase angles look good below 500 Hz. 

The source-solution SPL values are shown Fig. 5e and the source- 
solution location values are shown in Fig. 5f. Again, the variation of 
the source-A location solution below 164 Hz corresponds to the 
change in location of a low-frequency noise source from near the end 
of the potential core of the jet toward the nozzle as higher-frequency 


solutions are found. Consequently, we again associate SPL A with the 
jet source in the frequency range below 164 Hz. 

Again, the variation of the source amplitude of the solution for 
source B shown in Fig. 5e and the variation of the far-field SPL B 
autospectra shown in Fig. 5a have the same frequency dependence 
that the combustion modes have, in that they seem to be peaking in 
the 0- to 200-Hz band and the 200- to 400-Hz band. This 
characteristic behavior was reported in a study of PW4098 annular 
combustor modes by Miles [3], This suggests that the correlated 
noise source B is associated with the combustor source in the 0- to 
400-Hz frequency range. 

The source locations of the two point sources shown in Fig. 4f 
have an interesting symmetry that will be discussed in the next 
section. 

4. Nl con = 2743 rpm 

The highest power setting examined herein is the AH corr = 
2743-rpm operating condition. The autospectrum at 120 deg is 
shown in Fig. 6a. Below 100 Hz, the uncorrelated noise source is not 
needed to represent the measurements. The noise radiation pattern is 
well-represented again below 100 Hz by a single point source, as 
shown in Figs. 6b and 6c. In these figures, the magnitude and phase of 
the cross spectrum between the 120- and 130-deg microphones are 
shown. The results shown are for the three-parameter couplet-source 
model (A, N, and x A ) described in Appendix C. As shown in Fig. 6a, 
at frequencies above 100 Hz, the computer solution included an 
uncorrelated noise source. The coherence shown in Fig. 6d is above 
the coherence threshold y^ n = 0.0065 below 100 Hz. Below 100 Hz, 
the coherence is above 0.05 . Consequently, the cross-spectrum phase 
angles look good below 100 Hz. 

Above 200 Hz, the coherence is that of two random signals, except 
at select points. Consequently, the phase angle shown in Fig. 6c 
preceding 200 Hz is very random. The autospectrum solution shown 
in Fig. 6a looks good because the solution is based on mostly the 
incoherent level SPL^. Below 100 Hz, the uncorrelated noise level 
was too low to plot, as shown in Fig. 6e. As shown in Fig. 6f, below 
100 Hz, the decrease in the value of the location of source A with 
frequency shows that source A can again be inteipreted as being due 
to low-frequency sound radiation from the jet, as was the case with 
the .t A results shown in Figs. 4f and 5f. Again, this corresponds to the 
change in location of a low-frequency noise source from near the end 
of the potential core of the jet toward the nozzle as higher-frequency 
solutions are found. 

V. Source Location 

To provide further information on the source-location parameter, 
statistical distribution studies were conducted. In addition, for the 
five-parameter cases, scatter plots of x A and x B indexed by the 
frequency sequence number were made. These results are discussed 
in this section. 

A. /VI,,,,, = 582 rpm 

Histogram and nonparametric kernel-density estimates [43] of the 
distribution of location parameter x A are shown in Fig. 7. This plot 
also shows that the computer tends to find solutions in which 
x A = 10 ft at low frequencies and 20 < x A < 40 at higher 
frequencies. The nozzle in this geometry appears to be near 
x A = 10 ft. 

The greater phase-angle standard deviation in the frequency range 
above 400 Hz has increased the uncertainty of the location parameter 
for source A, x A . 

B. /VI,.,,,, = 1622 rpm 

A scatter plot of x A and x B identifying frequency with the index 
number is shown in Fig. 8a. The movement of source location x A 
from near the end of the jet potential core toward the nozzle as the 
frequency is increased is clearly shown. At low frequencies, the two 
major sources are the core noise and the jet noise. Above 400 Hz, the 
coherence becomes more irregular and the cross-spectrum phase 
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Location of source A, ft 

Fig. 7 Histogram and nonparametric kernel-density estimates of the 
distribution of location parameter x A , using the default bandwidth 
(heavier line) and half the default bandwidth (lighter line) based on a 
three-parameter model; A one-dimensional scatter plot (or rug plot) is 
shown at the bottom of the graph; ,Vl,.„ rr = 582 rpm. 


becomes more random. Consequently, the location results are more 
unreliable. They are shown in the scatter plot but are not very useful. 

Figures 8b and 8c show a histogram and nonparametric kernel- 
density estimates of the distribution of location parameter x A using 
the default bandwidth (heavier line) and half the default bandwidth 
(lighter line) for 0-152 and 164-937 Hz. A one-dimensional scatter 
plot (or rug plot) is shown at the bottom of the graphs. For the 
frequency range of 0-152 Hz shown in Fig. 8b, the distribution is flat 
with a concentration near 40 ft. For the frequency range of 164- 
937 Hz shown in Fig. 8c, the distribution has peaks between 1 1 and 
15 ft. 

Figures 8d and 8e show a histogram and nonparametric kernel- 
density estimates of the distribution of location parameter x B using 
the default bandwidth (heavier line) and half the default bandwidth 
(lighter line) for two frequency ranges: 0-152 and 164-937 Hz. A 
one-dimensional scatter plot (or rug plot) is shown at the bottom of 
the graphs. For the frequency range of 0-152 Hz shown in Fig. 8d, 
the distribution peaks between 1 1 and 15 ft. For the frequency range 
of 164-937 Hz shown in Fig. 8e, the distribution has one peak near 
—5 ft and another near 7 ft. 

C. /V 1 corr = 1900 rpm 

A scatter plot of x A and x B identifying frequency with the index 
number is shown in Fig. 9a. The movement of source location x A 
from near the end of the jet potential core toward the nozzle as the 
frequency is increased is clearly shown. At low frequencies, the two 
major sources are the core noise and the jet noise. Again, above 
400 Hz, the coherence becomes more irregular and cross-spectrum 




Fig. 8 Five-parameter model at 150 ft with )Vl c „ rr = 1622 rpm: a) scatter plot of x A by x B coded by frequency index for the whole region (164-937 Hz), 
b) histogram and nonparametric kernel-density estimates of the distribution of location parameter x A , (0-152 Hz), c) histogram and nonparametric 
kernel-density estimates of the distribution of location parameter x A (164-937 Hz), d) histogram and nonparametric kernel-density estimates of the 
distribution of location parameter x B (0-152 Hz), and e) histogram and nonparametric kernel-density estimates of the distribution of location parameter 
x B (164-937 Hz). 
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Fig. 9 Five-parameter model at /V 1 c „ rr = 1900 rpm: a) scatter plot of x A by x B coded by frequency index for the whole region (164-937 Hz), 
b) histogram and nonparametric kernel-density estimates of the distribution of location parameter x A (0-152 Hz), c) histogram and nonparametric 
kernel-density estimates of the distribution of location parameter x A (164-937 Hz), d) histogram and nonparametric kernel-density estimates of the 
distribution of location parameter x B (0-152 Hz), and e) histogram and nonparametric kernel-density estimates of the distribution of location parameter 
x B (164-937 Hz). 


phase becomes more random. Consequently, the location results are 
more unreliable. They are shown in the scatter plot but are not very 
useful. 

Figures 9b and 9c show a histogram and nonparametric kernel- 
density estimates of the distribution of location parameter x A using 
the default bandwidth (heavier line) and half the default bandwidth 
(lighter line) for two frequency ranges: 0-152 and 164-937 Hz. A 
one-dimensional scatter plot (or rug plot) is shown at the bottom of 
the graphs. For the frequency range of 0-152 Hz shown in Fig. 9b, 
the distribution is flat with a concentration near 30 ft. For the 
frequency range of 1 64-937 Hz shown in Fig. 9c, the distribution has 
one peak near 20 ft and another near 13 ft. 

Figures 9d and 9e show a histogram and nonparametric kernel- 
density estimates of the distribution of location parameter x B using 
the default bandwidth (heavier line) and half the default bandwidth 
(lighter line) for two frequency ranges: 0-152 and 164-937 Hz. A 
one-dimensional scatter plot (or rug plot) is shown at the bottom of 
the graphs. For the frequency range of 0-152 Hz shown in Fig. 9d, 
the distribution peaks near 15 ft. For the frequency range of 164- 
937 Hz shown in Fig. 9e, the distribution has one peak near —7 ft and 
another near 7 ft. 


D. /Vl t „ rr = 2743 rpm 

Histogram and nonparametric kernel-density estimates of the 
distribution of location parameter x A are shown in Fig. 10. This plot 
also shows that the computer tends to find solutions in which 


x A = 10 ft at low frequencies and 20 < x A < 40 at higher 
frequencies. The nozzle in this geometry appears to be near 
jf A = 10 ft. 

The phase-angle randomness in the frequency range above 100 Hz 
has increased the uncertainty of the location parameter for source A, 

x A . 

VI. Discussion 

The solution procedure provides location and sound-pressure- 
level values. The results just presented indicate that one can 
use the model solutions to separate the jet noise and the combustion 
noise source provided the coherence between the micro- 
phones is high enough to reduce the standard deviation of the 
phase angle. 

For the (Vl corr cases of 1622, 1900, and 2743 rpm, the procedure 
tracks the jet-mixing-noise source location as a function of 
frequency. These results are in good agreement with those obtained 
using linear phased microphone arrays to study the source 
distribution of a Rolls-Royce/SNECMA Olympus engine, as 
presented by Flynn and Kinns [ 27 1 (Figs. 6 and 7) and by Billingsley 
and Kinns [28] (Fig. 13). 

In addition, for the AH corr cases of 1622 and 1900 rpm, the 
procedure identifies the core-noise source at the nozzle exit. These 
results are in good agreement with those using the polar correlation 
technique of Fisher et al. [ 29 ] (Figs. 15 and 16), in which a source 
distribution of a Viper 601 is shown and in which a source 
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Location of source A, ft 

Fig. 10 Histogram and nonparametric kernel-density estimates of the 
distribution of location parameter x A , using the default bandwidth 
(heavier line) and half the default bandwidth (lighter line) based on a 
three-parameter model; a one-dimensional scatter plot (or rug plot) is 
shown at the bottom of the graph; AT,.„ rr = 2743 rpm. 


distribution of a RP21 1 is shown. In addition, they are similar to the 
results shown by Siller et al. [ 22 ] (Fig. 17), who used a phased- 
microphone array to study BR700 engine noise. 

These results show how important it is to start with no 
preconceived idea of what acoustic-source mode best represents the 
measured turbojet core noise. One should use a model that best fits 
the data and not one’s idea of the model. 

The results show that the uncorrelated noise source is important. 
Furthermore, the results show that the coherent core noise is greater 
than the coherent jet noise at frequencies less than 400 Hz at these 
angles. 

The use of four microphones separated by 10 deg in the method 
makes the procedure applicable to many sets of data taken in the 
normal study of engine noise at engine test stands. This makes this 
source-separation procedure easier to apply than other microphone- 
array methods that require special close microphone spacing and 
many more microphones [22,28,29,27] or the binaural source- 
location procedure of Kinns [3 1 ]. which requires close microphone 
spacing. In addition, the source-separation method might be 
applicable to the core-noise-source modeling portion of the Aircraft 
Noise Prediction Program (ANOPP) [44] . 


VII. Conclusions 

A method to separate correlated signals buried in turbofan engine 
core noise was presented. The method is based on finding acoustic 
model coefficients that enable a system of equations based on one or 
two correlated monopole noise sources and an uncorrelated 
monopole noise source to reproduce experimental data. 

At the lower engine-power settings, the nozzle exit is the strongest 
noise source. At the higher engine-power settings, the jet noise is 
observable, but the core noise is not observable because the 
coherence between the microphone becomes too low and the phase- 
angle measurements become meaningless. Above 164 Hz, at the 
higher engine-power setting, one has only an uncorrelated noise 
source. 

The effectiveness and the reliability of the method have 
been tested using Pratt & Whitney PW4098 far-field acoustic 
measurements. The agreement achieved between the experimental 


data and the acoustic model used confirm the effectiveness of the 
procedure. 


Appendix A: One-Sided Spectrum Functions 
for Acoustic Signals from a Point Source 

We shall assume a compact source region and assume that the 
source produces a wave that spreads spherically outward with no 
preferred direction. The wave equation in this case, as discussed by 
Morse and Ingard [45] (page 309), is 


i_a_( 2 ?p\ = ±Pp_ 

r 2 9rV dr) c 2 a dt 2 


(Al) 


We shall assume a simple source solution to the wave equation: 


P(t) = -exp (-ikR) sin(27r f 0 t) 
K 


where 

k = ~ c 2 =— -m/s 91 = 8314 J/kmol K (A2) 

c„ MW 


MW = 28.97 kg/kmol TF = 80 (A3) 

To get the autospectrum and cross spectrum of a point source, we use 
the following relationships from correlation and spectral analysis 
texts by Bendat and Piersol [4-6]. 

The autocorrelation function of x is 


R xx (r) = E[x k (t)x k (t + t)] (A4) 

The autocorrelation function for a sum of two processes y(t) = 
x { (t) + x 2 (t) is 

Ryy (t) = E\[Xl(t) + x 2 {t)][x l {t + t) + x 2 (t + t)]} 

= E XlX] (t) + R XiX2 (t) + R x , Xl (r) + R XlXl ( r) (A5) 


The two-sided auto- and cross-spectral density functions are 


S XX (J) = E[R xx (t)] = f R xx { r) exp(-;2jr/r) dr 

J— OO 

S X y(f) = F[R X y( r)] = f R xy ( r) exp(-;'2jr/r) dr 

J— OO 


(A6) 


and of course we have the inverse relationships 


/ OO 

s xx (f) exp(+y2rr/r) d / 

-OO 

r oo 

Rxy(r) = E~'[S xy (f)] = / S xy (f) expi+jlnfr) df 

J — OO 


(A7) 


The two-sided autospectrum function for a sum of two processes 
y(t) =x 1 (t) +x 2 (t) is 


Syy(f) = F[E{[x t (t) + x 2 (t)][xi(t + r) + x 2 (t + r)]}] 

= E[R XlXl (t) + R XiX 2 (t) + R X2 Xi ( r) + R X 2 j 2 ( r )] 

= S XlXl (/) + S XlXl (/) + S X2XI (/) + S X2X2 (f) ( A8) 


The autocorrelation function is the Fourier transform of the 
correlation of x{ t ) with itself and is related to the Fourier transform of 
x(t) by 


S xx (f) = F[R xx (t)] = f R xx ( r) exp(-/27r/r) dr = X(f)X*(f) 

J —OO 

(A9) 
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The cross-correlation function is the Fourier transform of the 
correlation of x(t) and y(t) and is related to the Fourier transform of 
x(t) and y(t) by 

S X y(f) = m*(T)] = f° R xy ( r) exp(— y'2jr/r) dr = X(f)Y*(f) 

J — OO 

(A10) 

The one-sided spectral density function is G x (f) = 2S x (f) 
0 < f < oo. 

The autocorrelation function of a sine wave x(t) = A cos(2jr/ 0 f) 
is 


A 2 

R X A?) = — cos(2jr/ 0 r) (All) 

The power spectral density of a sine wave x(t) = A cos(2jt f Q t) is 

G xx (f) = Y S (f~fo) (A12) 

The one-sided spectrum function for a sum of two processes y(t) = 
.xyff) + x 2 (t) is 

Gyy(f) = G XlXl (f) + G X1X2 (f) + G X2XI (/) + G X2X2 (/] (A13) 


Next, these relationships are applied to signals produced by acoustic 
waves from a point source and received at microphones i and j. 


A 

PaM ) = —exp (ikR a .) sm(2nf 0 t) 

Oj 

Paj(t) = — explifc^.) sin(27r/ o 0 

K a, 


(A 14) 


where R a . is the distance from the point source a to microphone i, and 
R a . is the distance from the point source a to microphone j. 

For the one-sided autospectral density and the one-sided cross- 
spectral density, we have 

G PatPa m=\^J(J~h) (A15) 


PajPu 


(f) = 


1 A 2 


2 R„.Rn 


„ jkiR., Vwr- 


S(f-fo) 


(A16) 


Appendix B: Correlated Noise Source A 

Consider the ith and y'th microphones and a single acoustic point 
source. As shown in Appendix A, the autospectra are 


1 A 2 

?! 


'Pa i Pa l (f)=2W S(f - f ° ) 


(Bl) 


G P ' lPa m =~ w - h ) 


(B2) 


and the cross spectrum is 


.(/) = ; 


1 A 2 


_ e -MR ai -R aj ) 8(f _ fo) 


(B3) 


autospectra are 


1 A 2 IN 2 

G„ ,.(f) = ^S(f-f 0 )+&^ r S(f-f 0 ) 

l K„ l K,, 


Pa< Pat ' 


1 A 2 


1 N 2 


(/) = -nr&if ~ /o) + - fo) 


and the cross spectrum is 


1 A 2 


' Pa, P a 




(Cl) 

(C2) 

(C3) 


Appendix D: Correlated Noise Source A 
and Correlated Noise Source B 

Consider the ith and y'th microphones, a point acoustic source A 
and a point acoustic source B. As shown in Appendix A, the 
autospectra are 


Gp.Pa, W = - / 0 ) + \~S(J - /„) 


(Dl) 


Gp.,p., 00 = - fo) + 2 g-Hf ~ fo) 


1 B 2 

3, 


(D2) 


and the cross spectrum is 
1 A 2 


Gp lP ,(f) = x 


-jk(R a .-R a .) 


P,P^’ 2R R 


S(f~fo ) 


- K b, K b, 


(D3) 


Appendix E: Correlated Noise Source A, Correlated 
Noise Source B, and Uncorrelated Source N 

Consider the ;th and y'th microphones, a single acoustic point 
source A, and a noise source N. As shown in Appendix A, the 
autospectra are 


G Pa,Pa, (/) - ^ + 


2 R 2 ^(/-/o) + ^5(/-/o) 


11V 2 

2Rl 


(El) 


1 A 2 


'p.,p., (f) = 2W S(f ~ /o) + \ w 8(f ~ fo) + ^ (f ~ /o) 


1 N 2 

2 Yl 


(E2) 


and the cross spectrum is 


G PiPj ( f ) = ^ e~ MR -‘ - fo) 


1 B 2 


_ e -MR bl -R bi ) s(f _ fo) 


(E3) 


Appendix C: Correlated Noise Source A 
and Uncorrelated Source N 

Consider the ith and y'th microphones, a single acoustic point 
source A, and a noise source N. As shown in Appendix A, the 


Appendix F: Cost Function 

The cost function described is for the third model, which is the 
five-parameter case [SPL^/), SPL A (/), SPL B (f), x A (f), and 
x B (f) |. Measured values of the autospectra splayy(f , t + 2,1 + 2), 
cross-spectra magnitude splayy(f,l + 2, m + 2), cross-spectra 
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Table FI Cost-function angle relationships 


Microphone identification 

#id. deg 

^model’ deg 

1 

100 

280 

2 

110 

290 

3 

120 

300 

4 

130 

310 


phase angle phay(f, £ + 2, m + 2), and coherence y 2 (f,£ + 
2, m + 2) are assumed to be available. The cost function is 
minimized at each frequency. 

The third model in the sequence assumes that the data are the result 
of a source triplet consisting of a coherent source at location x A (f ) 
with y A = 0.0 of strength SPL A (/) (identified with index 5), a 
coherent source at location x B (f) withy B = 0.0 of strength SPL fi (/) 
(identified with index 6), and an incoherent source of strength 
SPL N (f) at the location x N = 0.0 and y N = 0.0 (identified with 
index 7). 

The far-field microphone model angles used in the cost function 
are shown in Table Fl_ and indices 1 through 4 are used. The pressure 
parameters are transformed to pressures using 

Pa=A= r oPo 10 sp W 20 p B = B = r oPo 10 SPL ^ 20 
p N = N = r 0 p o 10 SPLw/2 ° 


The locations of the far-field microphones and sources are given by 


xg(l) = r 0 150 cos[# mmlel (f )tt/ 180] 
yg(i) = r 0 1 50 sin[0 model (f ) jr/ 1 80] 
xg(5) = r 0 x A yg{ 5) = 0.0 
yg(6) = 0.0 xg(l) = 0.0 


l = 1 , 2 , 3, 4 
1= 1,2,3. 4 
xg( 6) = r 0 x B 

yg( 7) = o.o 


(F2) 


The separation distance between the sources and the far-field 
microphones are given by 

r( 5, m) = V[*g(5) - xg(m)f + [yg(5) - yg(m)] 2 
m = 1,2,3, 4 

r( 6. m) = y/[xg( 6) - xg(m)f + [yg(6) - yg(m)] 2 

(F3) 

m = 1,2, 3, 4 

r( 7. m) = vW(7) - xg(m)] 2 + [yg(7) - yg(tn)f 
m = 1,2, 3, 4 


The autospectra at each microphone are calculated using 

tgy(l,l ) = tgu(£,l) + tgnu(i,l ) + tgnn(l,l) 

r 1 (F7) 

tsply(l, £) = 10.01og| 0 |^gy(f , (,)/ 

where l goes from 1 to 4. 

The autospectra error cost function suml is then 

95 0 4 

sum 1 = — ^[tsply(£,l) — splayy(f,l + 2, i 4- 2)] 2 (F8) 

® 1 = 1 

The cross-spectra magnitude and phase at each microphone are 
calculated using 

tgy(l , mi) = tgu(l, m) + tgnu(l, m) 

tsplyd , m) = 10.01og, 0 [fgy(f , (F9) 

tphay(l , m) = arctan{ ,3[igy(f , m)], ;H[tgy(f , wz)]} * 180/ 7 r 

where i goes from 1 to 3, and for each i value, m goes from l + 1 to 
4. 

The cross-spectra magnitude error cost function sum2 is then 

sum 2 = Yy sply ^' m) ~ s P la yy(f' (- + 2,m + 2)] 2 

(F10) 

where in doing the summation, l goes from 1 to 3, and for each t 
value, m goes from t + 1 to 4. 

In calculating the cross-spectrum phase-angle cost function, we 
first calculate a normalizing weighting factor from the measured 
coherence values: 

sumgam = ^ y 2 (f, l + 2, m + 2) (FI 1) 


where in doing the summation, l goes from 1 to 3, and for each l 
value, m goes from l + 1 to 4. 

Next, we calculate a phase-angle error based on the distance 
between points on a unit circle weighted by the measure coherence: 


angx(i, m) = cos [tphay{l, m) x 77 / 180] 
angy(l , m) = sin [tphay(l, m) x 77 / 180] (F12) 

angx 2(1, m) = cos [phay(f, l + 2, m + 2) x tt/180] 


The autospectra components are calculated as 

tgu(l , l) = 0.5p^/[r(5, l) 2 ] tnu(l , f) = 0.5p|/[r(6, l) 2 ] 

tgnn(t,€) = 0.5p 2 /[r{l,ty] (F4) 

where 1= 1,2, 3, 4. 

The cross-spectra phase is calculated using 

<p a (l, m) = — &[>-(5, l) — r( 5, m)] 

(F5) 

m) = —k[r( 6 , i) — r( 6 , m)\ 

where £ goes from 1 to 3 , and for each £ value, m goes from £ + 1 to 4 
and the wave number k is 2nf / c D . 

The cross-spectra phase components are 

tgu(l, m) = 0.5{exp[. j(j> a (l, m)]}p 2 A /[r( 5, l)r( 5, m )] 
tnu(£, m) = 0.5{exp[/^(f , m)\}p 2 B /[r(6, £)r( 6, m)\ 

where £ goes from 1 to 3, and for each £ value, m goes from £ + 1 to 
4. 


cingy2(l, m) = sin [phay(f, £ + 2, m + 2) x tt/180] 
angrad2 {£ , m) = [ angx(l , m) — angx 2(£, m)] 2 
+ [angy (£ , m) — angy2(l, m)] 2 
phawt{l,m) = y 2 (f,£ + 2, m 4- 2)angrad2(£, m) / sumgam 

(FI 3) 

The cross-spectra phase-angle-error cost function sum3wt is then 
sum 3wt = 1000 phawt(l, m) (FI 4) 

where in doing the summation, £ goes from 1 to 3, and for each £ 
value, m goes from £ + 1 to 4. 

The total cost function is 

sumf = suml + sum2 + sum3wt + conp2n50 

+ wt 1 + wt2 + wt3 (F15) 

where 
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conp2n50 - | _l ( 10 _ SPL^) 2 : if SPL W < 10 (F16) 

x5hi A = | 2 ° 0 0 0 i{ 164 x51owA = 1 1.0 (F17) 


x5hifi = I I°° n ° x51owB = -25.0 

[ 11.0 if / > 164 


I 01 

wt 1 = < ™ (x A - x5hi) 2 : if x A > x5hi 

[ ™ (x A — x51ow) 2 : if x A < x51ow 


(° : 

wt2 = < jjr (x B — x5hifi) 2 : if x B > x5hiB 

[ fr (x B — x51owS) 2 : if x B < x51ow5 


Wt 3 { (x B ~ x A ) 2 : if x A > x B 

The functions conp2n50, wt 1, and wt2 are used to keep the 
parameters within bounds. For / < 164, we have 1 1 < x A < 100 and 
—25 < x B < 100, and for / > 164, we have 1 1 < x A < 25 and 
—25 < x B < 1 1. The wt3 function is used to keep point source A 
separate from point source B. It removes the symmetry in the system 
model. The symmetry permits the interchange of the parameters 
associated with the coherent sources A and B without any change in 
the error. Rather than sorting the answers after the solution is 
obtained, the wt3 cost function forces the desired result. 
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